A General Phase Matching Condition for Quantum Searching Algorithm 
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A general consideration on the phase rotations in quantum searching algorithm is taken in this 
work. As four phase rotations on the initial state, the marked states, and the states orthogonal to 
them are taken account, we deduce a phase matching condition for a successful search. The optimal 
options for these phase are obtained consequently. 
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Quantum mechanical algorithms have recently become 
more and more popular in the field of computation sci- 
ence because they can speed up a computation when 
compared with classical algorithms. A famous example is 
the factorizing algorithm discovered by Shor|lJ . Another, 
which is what we intend to deal with in this work, is the 
quantum search algorithm developed by Grover 1| • If 
there is an unsorted database containing N items, and 
out of which only one marked item satisfies a given con- 
dition, then using Grover's algorithm one will find the 
object in 0(y/~N) quantum mechanical steps instead of 
O(N) classical steps. In many following researches it 
has been pointed out that the Grover algorithm is opti- 
mal. Benett et al.[[| showed that no quantum algorithm 
can search an object better than 0(yN) steps. Boyer et 
al.(^] gave tight bounds on Grover's algorithms. Zalka || 
has improved the tight bounds to show that Grover's al- 
gorithm is optimal and proposed a further improvement 
on the algorithm. Basically, Grover's algorithm consists 
of three successive unitary operations. They are (1) the 
Walsh-Hadamard transformation on the qubit |0) to cre- 
ate an initial state in the superposition of the basis states; 
(2) the 7r-angle rotation of the marked state; and (3) the 
inversion of the initial state. 

Although Grover [0 also proposed that the Walsh- 
Hadamard transformation can be replaced by almost any 
unitary transformation to create the initial state, in this 
work, however, we will focus on the phase rotations for 
the initial, the marked states, and the states orthogonal 
to them. Long et al. [8] have shown that when arbitrary 
rotations for the marked and the initial states are used 
instead of the inversions in Grover's original algorithm, 
the rotation phases must satisfy certain matching condi- 
tions. Galindo and Martm-Delgado Q recently have ad- 
dressed a more general viewpoint on the phase rotations. 
They gave four parameters in the Grover kernel since the 
rotations are operated on the four states including the 
initial, the marked, and the two states orthogonal to the 
former two. As the 7r-radian inversions for the initial and 
the marked states were considered, Galindo and Martin- 
Delgado concluded that the other two rotation phases 
must be equal. In this work, we will also study the four 
phase rotations in the Grover kernel but derive a gen- 
eral matching condition between the parameters without 



fixing any of them in advance. 

Suppose we have to search M objects out of N un- 
sorted elements, or in expressions: 



f( Wl ) = 1 ,4 = 1,2, ...,M. 
f( Wl ) = 0, i = M+l,M + 2, 



(1) 



,N. 



To begin the search, first construct a space spanned by 



the orthonormal set \w 



evenly distributed state as our initial state: 



1,2, ...,N, and give the 




(2) 



where 



w 



N 

E 



i 



y/N - M 



i=M+l 



and 



are unit and orthogonal 



Clearly, now both 

to each other. 

Next, let the Grover kernel G = — G2G1 operate on 
\s). Four phase parameters a, 0, 7, and S are given in 
the two Grover operators G\ and Gi respectively, such 

that in the space spanned by w ) and r ) , we have 



Gi 



a w 



■P r 



a 

(5 



g 2 = 7 |*><*i+<y(/-i*x«i) 



5 + (7 - S) M 



N 



(j-S) 



y/M(N-M) 

JV 



7-(7"^)f 



2 



Consequently, the Grover kernel then is given by 



G = - 



y/M(N-M) 



a(* + (7-*)%) /9(7-5)— - 



(3) 



Note that alternatively the parameters are written 



a = e Wl , (3 



= e^ 1 and 5 = e i<t>2 



1 , and <j) 2 are for the marked 



where the phases 6\, #2 
state, the state orthogonal to the marked state, the ini- 
tial state, and the state orthogonal to the initial state, 
respectively. 

Finally, repeat m times the operation, and wish the 

w) be greater 



probability of finding the marked element 
than some particular value. We then require, say, 



and then (w\gi) (gi\s) — 0(-^=), meaning that the re- 
quirement (4) will never be satisfied. Figure 1 shows an 
example for this case. To avoid (8), we in turn must have 
the matching condition 



a 7 
= 1 



or 9i — 62 — 4>\ 



(10) 



such that the eigenvector then becomes, in the normal- 
ized form 



l.9i> = 



13)8 



6)r 



/V2. 



We then have, under the matching condition (9), 



, , v , , x , M N-M^ . 1 
(w\9i)(9i\s) = (J- + ] /^^)/2^-. 



V = 



w 



G m \s) 



2 1 
> 2- 



(4) 



The amplitude (w G m \s) is deduced 



w 



G m \s) = + - l)Hsi> (5) 



where £ 12 and 1 171,2) denote the eigenvalues and the cor- 
responding eigenvectors of G, respectively. The detail 
expression of the eigenvalues and eigenvectors are given 
by 



£1,2 = e iAl ' 2 - ^TrG± ^(TrG) 2 - 4DetG, (6) 



1.91,2) oc 



where 



TrG 

DctG 
k 



fc±JV^/(TrG) 2 -4DctG 
2a(f-S)y/M(N-M) 
1 



M(a -/3)(j-S) + N(j(3 - aS) 
N ' 

N (7/3 - aS) - M(a + /3)( 7 - S). 



(7) 



(8) 



The eigenvector \gi) should be discussed in detail be- 
cause, as shown in (5), it determines whether the require- 
ment (4) can be satisfied. As M « N and N » 1, the 
eigenvector \gi) can be asymptotically expressed by 



Now, since 



- 1| = w 



1AA 



mAA 



sin ■ 



where AA = Ai — A 2 , the probability then is 



w 



G m \s) 



1AA 



+ (e* mAX -l)(w\g 1 )(g 1 \s} 



sin ■ 



(11) 



Figure 2 shows the p vs. m diagram under the condition 
(9). It shows that p can approach unity as |toAA| = 
(2j — l)n, j = 1, 2, .... Taking first n , we have 



m 



AA 

7T fW 

2 V ~N 

TT [~M 

2 V iV 



\8 



( a -0)(rr-6) 
1 

2cos(<j) 1 - 4> 2 ) 



(12) 



Eventually, the searching time step will be minimum 
when 



01-6-2 = TT, 

and the corresponding time step is 



l.9i) cc 



aS— 0j_ / N 
a(j-S) V M 
1 



II! 



(9) 



min(m) = 



TT M 

4 V ~N' 
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In figure 3, the p vs. m diagram is shown under an opti- 
mal choice of the parameters. 

To summarize, as the four parameters a, (3, 7, and 8 
are taken into consideration in a generalized Grover's ker- 
nel G, the phase matching condition for these parameters 
has been deduced. It is found that the marked elements 
will never be searched unless aS — (3j — 0. That is, the 
eigenvector \gi) of the Grover kernel G should not asymp- 
totically coincide the marked state. The optimal option 
for the relation between these parameters, however, is 

- = J = 1 or 0, - 9 2 =4> x - <P 2 = tt (13) 



The choice of parameters taken in the original Grover 
operator is only the simplest one, in which a = 7 = — 1 
and [3 = 8 = 1 were used. Long et al. [|| have treated 
the phase rotations for the initial and the marked states 
alone, and the phase matching condition is exactly the 
case of a — 7 and (3 = 8=1. Galindo and Martin- 
DelgadoQ, on the contrary, have discussed the phase 
rotations for the states orthogonal to the initial and the 
marked states, and have addressed the condition a = 7 = 
— 1 and (3 = 8. 



[1] P. W. Shor, in Proceedings of 35" 1 Annual Symposium on 
Foundation of Computer Science, (IEEE ,Los Alamitos, 
CA, 1994). 

[2] L. K. Grover, in Proceedings of 28 t/l Annual ACM Sym- 
posium on the Theory of Computation, (ACM Press, New 
York, 1996), p. 212. 

[3] L. K. Grover, Phys. Rev. Lett. 79, 325(1997). 

[4] C.H. Bennett, E. Bernstein, G. Brassard, and U. Vazirani, 
SIAM J. Comput. 26,1510(1997). 

[5] M. Boyer, G. Brassard, P. H0yer, and A. Tapp, Fortsch. 
Phys.-Prog. Phys., 46, 493(1998). 

[6] C. Zalka, Phys. Rev. A, 60, 2746(1999). 

[7] L. K. Grover, Phys. Rev. Lett. 80, 4329(1998). 

[8] G. L. Long, Y. S. Li, W. L. Zhang, L. Niu. Phys. Lett. A 
262, 27 (1999). 



[9] A. Galindo and M. A. Martm-Delgado, Phys. Rev. A 62, 
062303 (2000). 



Figure Captions 

Fig.l Probability p as a function of the time step m 
for N = 1000, M = 10, a = 7 = e l7r , (3 = e l % , 8 = ie a . 

Fig. 2 Probability p as a function of the time step m 
for N = 1000, M = 10, a = e ll - r * , 7 = e t7r , (3 = e ll - 67T , 
8 = e 10 - 9 *. 

Fig. 3 Probability p as a function of the time step m 
for N = 1000, M = 10, a = e ilJ * , 7 = e* 1 - 9 ^, (3 = e i0 - 77T , 
8 = e 10 - 9 *. 



